Non-Gaussianity in Two-Field Inflation 
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We derive semi-analytic formulae for the local bispectrum and trispectrum in general two-field 
inflation and provide a simple geometric recipe for building observationally allowed models with 
observable non-Gaussianity. We use the SN formalism and the transfer function formalism to express 
the bispectrum almost entirely in terms of model-independent physical quantities. Similarly, we 
calculate the trispectrum and show that the trispectrum parameter tnl can be expressed entirely in 
terms of spectral observables, which provides a new consistency relation unique to two-field inflation. 
We show that in order to generate observably large non-Gaussianity during inflation, the sourcing 
of curvature modes by isocurvature modes must be extremely sensitive to the initial conditions, 
and that the amount of sourcing must be moderate in order to avoid excessive fine-tuning. Under 
some minimal assumptions, we argue that the first condition is satisfied only when neighboring 
trajectories through the two-dimensional field space diverge during inflation. Geometrically, this 
means that the inflaton must roll along a ridge in the potential V for some time during inflation and 
that its trajectory must turn slightly (but not too sharply) in field space. Therefore, it follows that 
two-field scenarios with attractor solutions necessarily produce small non-Gaussianity. This explains 
why it has been so difficult to achieve large non-Gaussianity in two-field inflation, and why it has 
only been achieved in a narrow class of models like hybrid inflation and certain product potentials 
where the potential and/or the initial conditions are fine-tuned. Some of our conclusions generalize 
qualitatively to general multi-field inflation. 



I. INTRODUCTION 



Cosmological inflation [lHa] is widely thought to be re- 
sponsible for producing the density perturbations that 
initiated the formation of large-scale structure. Dur- 
ing such an inflationary expansion, quantum fluctuations 
would have been stretched outside the causal horizon and 
then frozen in as classical perturbations. These primor- 
dial perturbations would later be gravitationally ampli- 
fied over time into the cosmological large-scale structure 
that we observe today [6Hll|. 

Pinning down the specific nature of inflation or what- 
ever physics seeded the primordial density fluctuations is 
one of the greatest open problems in cosmology. The sim- 
plest models of inflation are driven by a single scalar field 
producing fluctuations that are adiabatic, nearly scale- 
invariant, and nearly Gaussian, but these assumptions 
need to be tested. Whether the primordial fluctuations 
were indeed adiabatic and near scale-invariant can be de- 
termined by measuring the power spectra of fluctuations; 
the upper limit on the isocurvature spectrum constrains 
non-adiabaticity, while the slope of the scalar (curva- 
ture) power spectrum constrains the deviation from scale- 
invariance. Similarly, whether the primordial fluctua- 
tions obey Gaussian statistics can be tested by measuring 
reduced n-point correlation functions, where n > 3, since 
for Gaussian fluctuations, these higher-point functions 
all vanish and only the two-point function (the power 
spectrum) is non-zero. Any deviations from adiabatic- 
ity, scale-invariance, and Gaussianity would signal some 
non-minimal modifications to the simplest scenarios, and 
hence would provide exciting insight into ultra high en- 
ergy physics. 



Of these observational measures, non-Gaussianity has 
the potential to be the most discriminating probe, given 
all the information contained in higher-point statistics. 
This is particularly valuable given how challenging it has 
been to discriminate among the myriad different infla- 
tionary models. 

The two lowest order non-Gaussian measures are the 
bispectrum and the trispectrum. Just like the power 
spectrum V-r represents the two-point function of the 
comoving curvature perturbation 1Z in Fourier space, the 
bispectrum Btz represents the three-point function and 
the trispectrum T-r represents the four-point function: 

(ft(ki)ft(ka)) = (27r) 3 5 (j2 k ^J 7Mki,k 2 ), (1) 



<ft(ki)ft(k 2 )ft(ka)) - (2 7 r) 3 5(^k i j Bfc(ki,ka,k3 



vi=l 



(2) 



(7e(k 1 )7e(k 2 )^(k 3 )^(k 4 )) = 

= (2ir) 3 6 (j^^ xr TC (k 1 ,k 2 ,k 3 ,k 4 ). (3) 

The 5-functions in equations CD)-© reflect the fact that 
the statistical properties are translationally invariant in 
real space, which makes the above three correlations van- 
ish except if all k- vectors add up to zero — i.e., are the 
negative of one another for the power spectrum, form the 
three sides of a triangle for the bispectrum, and form the 
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sides of a (perhaps non-flat) quadrangle for the trispec- 
trum. Since the statistical properties are also rotationally 
invariant, the power spectrum depends only on the length 
and not the direction of its vector argument, and the bis- 
pectrum depends only on the lengths of the three triangle 
sides, so that they can be written simply as Vn(k) and 
Bu(kx, &2) fe), respectively. 

For the bispectrum, it is standard in the literature to 
define a dimensionless quantity /atl(A;i, fe, fe) to repre- 
sent the bispectrum by dividing by appropriate powers 
of the power spectrum |12j| 1 : 



-fNL(ki,k 2 ,k 3 ) 
5 



■Pn(h)Vn(k2)+ ,V " U ' 



permutations 



(4) 



Although /jvi, can in principle depend on the triangle 
shape of the k-vectors in very complicated ways, it has 
been shown that in practice, essentially all models pro- 
duce an /tvl that is well approximated by one of merely 
a handful of particular functions of triangle shape, with 
names like "local", "equilateral", "warm", and "flat" 
[HIE]- For example, the local function peaks around tri- 
angles that are degenerate (with one angle close to zero, 
like for k 3 -C &i ~ fca); while the equilateral function 
peaks around triangles that are equilateral (ki = hi = 
ks). Bispectra dominated by different triangle shapes 
correspond to different inflationary scenarios and differ- 
ent physics. In particular, for multi-field inflation, bar- 
ring non-canonical kinetic terms or higher-order deriva- 
tive terms in the Lagran gian , the dominant type of bis- 
pectra is of the local form Local non-Gaussianity 
arises from the non-linear evolution of density perturba- 
tions once the field fluctuations are stretched beyond the 
causal horizon. Seven-year data from WMAP constrains 
non-Gaussianity of the local form to [l5[ 



10 < Inl < 74 (95% C.L. 



(5) 



and a perfect CMB measurement has the potential to 
detect a bispectrum as low as |/jvx| ~ 3 [161 ]. 

Similarly, the dominant form of trispectra for standard 
multi-field inflation is also of the local form, and it can 
be characterized by two dimensionless non-linear param- 
eters, tnl and gNL- Five-year data from WMAP con- 
strains these two parameters to [TtJ 

-0.6 x 10 4 < t nl < 3.3 x 10 4 (95% C.L.), (6) 
-7.4 x 10 5 < g NL < 8.2 x 10 5 (95% C.L.). (7) 

Interestingly, for single-field inflation, the non-linear 
parameters representing the bispectrum [l2|, [l8l - l25l | and 



1 The factor of —2 arises from the fact that /jvl was originally 
defined with respect to the metric perturbation. After inflation 
ends and during the matter-dominated era, 2$ = -%7L. 



trispectrum [26H28fl are all of order the slow-roll param- 
eters (i.e., at the percent level) and will not be acces- 
sible to CMB experiments. However, if inflation is de- 
scribed by some non-minimal modification, such as multi- 
ple fields or higher derivative operators in the inflationary 
Lagrangian, then non-Gaussianity might be observable in 
the near future. Indeed, there have been many attempts 
to calculate the level of non-Gaussianity in general multi- 
field models (e.g.. Jg^M |3^-|io[ |6^). as well as in two- 
field models (e.g., (4l] - l49l |). However, it has been very dif- 
ficult to find models that produce large non-Gaussianity, 
though some exceptions have been found such as in the 
curvaton model |5CH55j|. hybrid and multi-brid inflation 
(e.g., j56l - [64| ). and in certain modulated and tachyonic 
(p)re-heating scenarios (e.g., [65l - l7l| ). Moreover, it has 
not been wholly clear why it is so difficult to produce 
large non- Gaussianity in such models. Though some au- 
thors [H HO, [M lH, |H EH have found spikes in non- 
Gaussianity whenever the inflaton trajectory changes di- 
rection sharply, these spikes in non-Gaussianity are tran- 
sient and die away before the end of inflation. That 
makes a comprehensive study of non-Gaussianity gen- 
eration timely, to understand any circumstances under 
which observably large non-Gaussianity arises in such 
models. 

In this paper, we calculate the bispectrum and trispec- 
trum in general inflation models with standard kinetic 
terms, focusing on the important case of two-field infla- 
tion. We provide conditions for large non-Gaussianity 
and a unified answer to the mystery of why it has been 
so hard to produce large non-Gaussianity in two-field in- 
flationary models. The rest of this paper is organized 
as follows. In Section [Hi we present the background 
equations of motion for the fields and discuss the field 
vector kinematics. Section Mil presents the equations of 
motion for the field perturbations and some necessary re- 
sults for the power spectra. In Section IIV1 we describe 
the 8N formalism, which we use to calculate the bispec- 
trum. In tandem, we discuss the necessary conditions for 
large non-Gaussianity. Finally, we tackle the trispectrum 
in Section [Vj We summarize our conclusions in Section 



II. BACKGROUND FIELD EQUATION & 
KINEMATICS 



In this section, we review the background equations 
of motion and discuss the kinematics of the background 
fields. This discussion will help us calculate the pri- 
mordial bispectrum and trispectrum in two-field infla- 
tion and understand what features are necessary for non- 
Gaussianity to be observably large. 

We consider general two-field inflation where the non- 
gravitational part of the action is of the form 



i d^df_ 

2 9 l3 dx»dx v 



-gd A x, (8) 
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where V(fa, fa) is a completely arbitrary potential of the 
two fields, is the spacetime metric, and &y reflects the 
fact that we assume the kinetic terms are canonical. 

In this paper, we adopt similar notation to [z2|: bold- 
face for vectors, the symbol T to denote the transpose 
of a vector, and standard vector product notation. Gra- 
dients and partial derivatives represent derivatives with 
respect to the fields, except where explicitly indicated 
otherwise. We set the reduced Planck mass, fh = 8wG, 
equal to unity, so that all fields are measured in units of 
the reduced Planck mass. To simplify the equations of 
motion and connect them more directly with observables, 
we use the number of e-folds, N, as our time variable. N 
is defined through the relation 



dN = Hdt, 



(9) 



where t is the comoving time and H is the Hubble pa- 
rameter. We denote derivatives with respect to N using 
the notation 



_d_ 

dN' 



(10) 



Using N as the time variable, we showed in |72| that 
the background equation of motion for the fields can be 
written as 



V 



<p' = -VlnF. 



(3-6) 

The parameter e is defined as 

e=-(lniJ)' = ±0' ■</>', 
and r] is the field acceleration, defined as 
V = • 



(11) 



(12) 



(13) 



In (72| , we also explained how the two quantities <f>' 
and r) represent the kinematics of the background fields. 
If we view the fields as coordinates on the field manifold, 
then <f> represents the field velocity, and 



v=W\ 



(14) 



represents the field speed. Similarly, rj is the field accel- 
eration. 

The velocity vector, <fi' , is also useful because it can be 
used to define a kinematical basis [73- 75]. In this basis, 
the basis vector eii points along the field trajectory, while 
the basis vector e±_ points perpendicularly to the field 
trajectory, in the direction that makes the scalar product 
e± ■ 77 positive. To denote the components of a vector and 
a matrix in this basis, we use the short-hand notation 



X 



'II 



and 



M||_l = ef Meu 



etc. 



(15) 



(16) 



The kinematical basis is useful for several reasons. First, 
the field perturbations naturally decompose into com- 
ponents parallel and perpendicular to the field trajec- 
tory, and the former represent bona fide density pertur- 
bations, while the latter do not. This decomposition of 
the field perturbations is helpful in finding expressions for 
the power spectra. Second, it allows us to consider sep- 
arate aspects of the background field kinematics, which 
in [72j, we encapsulated in a set of three quantities. The 
first quantity is the field speed, v. The second and third 
quantities arise from decomposing the field acceleration 
into components parallel and perpendicular to the field 
velocity. In particular, the quantity — represents the log- 
arithmic rate of change in the field speed (the speed-up 
rate), while the quantity ^ represents the rate at which 
the field trajectory changes direction (the turn rate) [72j. 

This distinction between the speed-up rate and the 
turn rate is important for two reasons. First, the 
turn rate represents uniquely multi-field behavior (as the 
background trajectory cannot turn in single-field infla- 
tion), whereas the speed-up rate represents single-ficld- 
like behavior. Second, the speed-up and turn rates have 
very different effects on the evolution of the field per- 
turbations and hence on the power spectra. Indeed, the 
features in the power spectra depend not only on the ab- 
solute sizes of the two rates but also on their relative sizes 
to each other; in particular, the ratio of the turn rate to 
the speed-up rate is an indicator of the relative impact 
of multi-field effects. So disentangling the two quantities 
allows for a better understanding of the power spectra 
and all the ways that the spectra can be made consistent 
with observations. 

To fully take advantage of this distinction between the 
speed-up and turn rates, we redefined the standard slow- 
roll approximation in 72], splitting it into two different 
approximations that can be invoked either separately or 
together. As background, the standard slow-roll approx- 
imation is typically expressed as 



e«i|VlnV| 2 «l, 



and 



didjV 



V 



< 1. 



(17) 



(18) 



However, as argued in [72j, the latter condition lumps to- 
gether and simultaneously forces the speed-up rate, the 
turn rate, and a quantity called the entropy mass to be 
small. So instead, we redefined the slow-roll approxima- 
tion to mean that the field speed is small, 



e=i, 2 «l ! 



and is slowly changing, 



3L<i. 



(19) 



(20) 
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In other words, the above slow-roll approximation rep- 
resents the minimum conditions necessary to guaran- 
tee quasi-exponential inflationary expansion and corre- 
sponds to limits on single-field-like behavior. As for the 
turn rate, we endowed it with its own separate approxi- 
mation, the slow-turn approximation, which applies when 
the turn rate satisfies 



V-i- 



< 1. 



(21) 



The slow-turn limit corresponds to limits on multi-field 
behavior. Finally, this alternative framework allows the 
lowest order entropy mass to take on any arbitrary value. 
(See [z3l for further discussion of these points.) 

When the background field vector is both slowly rolling 
and slowly turning, we call the combined slow-roll and 
slow-turn limits (which is equivalent to the conventional 
slow-roll limit minus the constraint on the entropy mass) 
the SRST limit for brevity. In the combined limit, the 
evolution equation for the fields can be approximated by 



4>' 



-VlnV. 



(22) 



Also in this combined limit, the speed-up rate and the 
turn rate can be approximated by 



m 

V 



-Mi 



v 



-Mi 



ll-L: 



(23) 



respectively, where we define the mass matrix, M, as the 
Hessian of In V, i.e., 



M = VV J mV. 



(24) 



Being a symmetric 2x2 matrix, M is characterized by 
three independent coefficients. In the kinematical ba- 
sis, these three coefficients are M^, M^ ± , and M±±, 
the third of which is the lowest order entropy mass. 2 In 
other words, in the kinematical basis and under the SRST 
limit, we can interpret the mass matrix as follows: 



M = 



Mini Mi 



Mi 



ll-L = 

M xx 



-speed-up rate —turn rate 
—turn rate entropy mass 



(25) 



where the speed-up rate and turn rate alone determine 
the background kinematics. However, all three quantities 

— the speed-up rate, the turn rate, and the entropy mass 

- affect how the perturbations evolve, as described in 
the next section. 



III. PERTURBATIONS, TRANSFER 
FUNCTIONS, AND POWER SPECTRA 

In this section, we summarize the general results for 
the evolution of perturbations and for the power spec- 
tra, as these expressions will enable us to calculate the 



bispectrum and trispectrum in two-field inflation and to 
express the results in terms of other spectral observables. 

In [72j], we derived the following general equation of 
motion for the field perturbations in Fourier space: 



(3- 



1 D 2 d(p D64> 



dN 



0IN 2 



a 2 V 



S4> 



M 



T 1 



(3- e ) 



S4>, (26) 



where 8<p represents the field perturbation in the 
flat gauge, which coincides with the gauge-invariant 
Mukhanov- Sasaki variable [7(1 |z3|. For modes in the 
super-horizon limit (k <C aH), we showed that when the 
background fields are in the SRST limit, the above ex- 
pression reduces to [721 ] 



DScf) 

dN 



-M5<f). 



(27) 



Now we switch to working in the kinematical basis, 
where the modes decompose into adiabatic modes, 5(j)\\ , 
and entropy modes, S(f>±. The former represent den- 
sity perturbations, while the latter represent relative field 
perturbations that leave the overall density unperturbed. 
In this basis, the super-horizon equations of motion for 
the two mode types are 



>J X * -Mx±6<l> ± , 



(28) 



where the first equation is exact and the second is valid to 
lowest order in the slow-turn limit. (Full expressions are 
given in (72[.) In the SRST limit, the evolution of modes 
is determined by the three unique coefficients of the mass 
matrix. The evolution of adiabatic modes is controlled 
by « -M|||| and by ^ w -M\\ ± . The third unique 
coefficient of the mass matrix, M±±, alone determines 
the relative damping or growth of entropy modes. We call 
Mj__l the lowest order entropy mass (or just the entropy 
mass) because it approximates the effective mass in the 
full second-order differential equation of motion for the 
entropy modes 72]. In addition to being viewed as an 
effective mass, M±± can also be viewed as a measure of 
the curvature of the potential along the ej_ or entropic 
direction. When the curvature of the potential along the 
entropic direction is positive, the entropy modes decay; 
when the curvature is negative, the entropy modes grow. 

Directly related to these two modes are the curvature 
and isocurvature modes, the two quantities whose power 
spectra are typically computed when considering the two- 
field power spectra. During inflation, the curvature and 
isocurvature modes are simply related to the adiabatic 
and entropy modes, respectively, by a factor of i [78j . 
That is, the curvature modes are given by 



2 We refer to M±± as the entropy mass, even though we con- 
structed it to be dimensionless. 



K 



'II 



(29) 
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and the isocurvature modes by 



S = 



(30) 



The super-horizon evolution of curvature and isocur- 
vature modes can be determined from the equations 
of motion for the adiabatic and entropy modes. We 
parametrize the solutions through the transfer matrix 
formalism WjlWsk: 



1 T nS \ 
T ss \ 



(31) 



where the transfer functions can be written as 

Tns(N*,N) = / a(N)T ss {N*,N)dN, 



T ss (N*,N) = eI», M d *. 



(32) 



The subscript * means that the quantity is to be eval- 
uated when the corresponding modes exit the horizon. 
The transfer function T$s therefore represents how much 
the isocurvature modes have decayed (or grown) after ex- 
iting the horizon. The transfer function T-rs represents 
the total sourcing of curvature modes by isocurvature 
modes; that is, it represents the importance of the multi- 
field effects. In [731, we found that 



2^ 



(33) 



exactly, which tells us that the curvature modes are only 
sourced by the isocurvature modes when the field trajec- 
tory changes direction. However, the isocurvature mass, 
/3, must be approximated or computed numerically. To 
lowest order in the SRST limit, 



a = -2M||_l, 

P = M U -M ±± . 



(34) 



From a geometrical perspective, equation AMI) for f3 
shows that how fast the isocurvature modes evolve de- 
pends on the difference between the curvatures of the po- 
tential along the entropic and adiabatic directions. From 
a kinematical perspective, the isocurvature modes will 
grow if the entropy modes grow faster than the field vec- 
tor picks up speed. This means that isocurvature modes 
tend to grow in two types of scenarios: when M±± is 
large and negative and when e decreases quickly. Other- 
wise, when the entropy modes do not grow faster than 
the field vector picks up speed, the isocurvature modes 
decay. 

Now we present expressions for the power spectra and 
their associated observables. The power spectrum of a 
quantity X is essentially the variance of its Fourier trans- 
form: 



(A'(k 1 )A'(k 2 )) = (2TT) 3 S(k 1 +k 2 )V x (k 1 ). 



(35) 



Using the above results, the curvature, cross, and isocur- 
vature spectra at the end of inflation can be written to 
lowest order as [7al 



Cns 



a* 

2^ 

H_ 

2tt 

H_ 

2n 



1 

2Z 



(1 



— I ~ I - — T-jisTss, 



2e* 



2e» 



L SS' 



(36) 



where it is implied that the transfer functions are eval- 
uated at the end of inflation. The associated curvature 



spectral index is [72 1 

d\nP n 



1 



dN 



[n T + 2elMe N ] r , (37) 



where n s is the standard scalar spectral index that is 
constrained by observations, ut = — 2e» is the tensor 
spectral index, and the unit vector e^v points in the di- 
rection of the gradient of N. In the kinematical basis, 
e^r takes the form [j^] 

e N = cos A N eji + sin A N e*j_, (38) 

where Ajv is the correlation angle, which satisfies 

tanAjv = 2W (39) 

The correlation angle can be given in terms of the dimen- 
sionlcss curvature-isocurvature correlation, rc , which we 
define as [721 



rc 



C. 



ns 



= sin A 



N, 



in analogy to the tensor-to-scalar ratio, 

V T 9 A 

r T = — = 16e*cos A Nl 
l J il 



(40) 



(41) 



where Vt is the tensor spectrum of gravitational waves. 
Similarly, another ratio of spectra, the isocurvature frac- 
tion, 



f. =^ = _ 



T. 



^ = cos 2 A w Tj 5 , (42) 



ns 



gives the relative size of T$s to T n $- 

These results provide the starting point for our com- 
putation of the bispectrum and trispectrum for two-field 
inflation. 



IV. THE BISPECTRUM 

In this section, we calculate the bispectrum for gen- 
eral two-field inflation using the 5N formalism and the 
transfer function formalism. We also develop a set of 
conditions encapsulating when the bispectrum is large 
enough to be detected. 
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A. Calculation of /jvl Using the SN Formalism 

We focus on bispectrum configurations of the local or 
squeezed type {e.g., k 3 <C fei ~ £2), which is the dom- 
inant type present during standard multi-field inflation. 
As we showed earlier in equation (j4]) , the bispectrum can 
be expressed in terms of the dimensionless non-linear pa- 
rameter Jnl [13] • 3 From here on, whenever Jnl appears 
in this paper, it represents the local form, so we drop the 
superscript local. 

Convenientl y, f /vT, can be written in terms of the SN 
formalism [8lMS3j] , where N again represents the num- 
ber of e-folds of inflation. Under the SN formalism, 
1Z = VJV • S(j> [13, where 5<fi is measured in the 
flat gauge and where it is implied that the gradient is 
with respect to the fields at horizon exit. (For brevity, 
we drop the subscript * on V in this section, but restore 
it in later sections whenever there might be some poten- 
tial ambiguity.) Using this result, correlators of 1Z can be 
written in terms of gradients of N. In particular, it has 
been shown [HI that the local form of /jvx can written 
as 

V T NVV T NVN 



6 f - 

'TJNL — 

5 



IVJV!- 



(44) 



We use equation (|44|) to find an expression for Jnl 
in terms of the transfer function formalism. We start 
by first finding a semi-analytic formula for VAT in two- 
field inflation. By comparing equation (|36|) to the lowest 
order result for the curvature power spectrum in multi- 
field inflation I82I1. 



2^ 



|VAf 



we obtain 



|VJV| = 



T 2 



2c, 



(45) 



(46) 



where again it is implied that T-jis is evaluated at the 
end of inflation. Combining equation (|46[) with the fact 
that 



4> ■ VN = 1, 



(47) 



we conclude that V7V takes the following form in the 
kinematical basis: 



VN 



/2el V 



(48) 



The above equation implies that we can also write ViV 
as 



VN 



1 



T 2 



2e» 



/2e* cos A 



N 



(49) 



where e^r is the unit vector in the direction of ViV and 
is given by equation (155)) . 

Next, we re- write equation (|4"4"|) for Jnl as 



6, 



e T N VV T Ne 



N 



ViV 



(50) 



Since e« ■ ejy = 1, it follows that Ve^v • &n = 0, and 
hence 



W'iVe 



N 



V|ViV| 



(51) 



Taking this result, dividing through by |VJV| and using 
equations (f3T)f and (f4*5)) . we find 



|ViV] 



2e* 



■sinAjvcosAjvVr TC<s . (52) 



In the SRST limit, using equations flT2j| and (|22|). it holds 
that 



Ve 



-M<f>. 



(53) 



Substituting this result into equation (1521) and dividing 
through by another factor of |ViV|, we find that 



VV 1 Ne N 
\VN\ 2 



= cos An x 



(54) 



M*e| + sin A^v cos A N ^/2cZ VT ns 



To complete our calculation of /nl, we need to 
contract equation (|54|) with the unit vector ejy. We 
break this calculation into two parts, based on the fact 
that ejy = cosAjve| + sinA^e^. First, we contract 

cos Ajv(e|) T with equation (|M)l . Using ^ = $ ■ V and 
the norm of equation (|49|) , we can write 



cosA^v (e|) T VV T iVe 



N 



|ViVP = cos 2 A N -j-^ In |ViV|. (55) 



From equations (l37l) and (|45|) . it follows that 
cos A N (e|) T VV T iVe 



iN 



1 



IVJVI 



cos An (n-!z — nx). (56) 



Second, we calculate sin An ( e l) T contracted with equa- 
tion (fM)) , which yields 



sinAjv {e*j_) T VV T N e N 



IVJVI 



sin A 7v cos A^x 



(57) 



A I 



II -L 



sin An cos A n v / 2eT • VT^ 



Combining equations (|56[) and (|57p . /nl can be written 
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3 The non-linear parameter /jyi was originally introduced to rep- 
resent the degree of non-Gaussianity in the metric perturbation 



(43) 



where <J><3 is Gaussian and <3? is not. Here, <I> is the metric pertur- 
bation in the Newtonian gauge, which equals the gauge-invariant 
Bardecn variable. 



Jnl = ~ (niz ~ tit) cos 2 A at + (^M*\± + sin An cos ANV2eZe* ± ■ VT^j sin A N cos A N . 



(58) 



Equation (|58l) depends on sines and/or cosines times 
the curvature (scalar) spectral index, the tensor spec- 
tral index, the turn rate 2i. ~ _ Mjn at horizon exit, 
and \/2e*e^ • VT^s. Observational constraints force the 
magnitudes of n-R. and tit to be much less than unity, 
and the turn rate at horizon exit must be at least some- 
what less than unity to avoid violating scale-invariance 
and causing a complete breakdown of the SRST approxi- 
mation at horizon exit. Therefore, the magnitude of Jnl 
cannot exceed unity unless 



I sin 2 A N cos 2 A N (e* ± ■ VT ns )\ > 



1 



(59) 



Physically speaking, for equation J5 
quires that two conditions be met: 



to be satisfied re- 



1 . The total amount of sourcing of curvature modes by 
isocurvature modes (Tus) must be extremely sen- 
sitive to a change in the initial conditions perpen- 
dicular to the inflaton trajectory. In other words, 
two neighboring trajectories must experience dra- 
matically different amounts of sourcing. 

2. The total amount of sourcing must be non-zero 
(i.e., sinA^ ^ 0). Usually, the amount of 
sourcing must also be moderate, to avoid having 
sin 2 A N cos A 2 N <gC 1. 



The first condition, that T-jis be extremely sensitive to 
the initial conditions, makes sense on an intuitive level. 
In order to produce large non-Gaussianity, perturbations 
off the classical trajectory must move the inflaton onto 
neighboring trajectories that experience very different 
dynamics for the field perturbations, hence producing a 
large degree of skew in the primordial fluctuations. Meet- 
ing this condition that Tus is very sensitive to the initial 
conditions can be achieved in a couple of different ways. 
Since equations (|52")l and show that Tus depends on 
an integral of the turn rate times the relative amplitude of 
isocurvature modes (Tss), neighboring trajectories need 
to have very different turn rate profiles, Tss profiles, or 
both. 

The second condition reflects the fact that in the limit 
of no sourcing — which corresponds to single-field be- 
havior — the bound in equation (|59l) can never be sat- 
isfied. Usually, the sourcing must also be moderate, 




FIG. 1: The trigonometric factor sin 2 Ajv cos 2 Ajv as a func- 
tion of the total amount of mode sourcing, Tus- For reference, 
when Tizs = 1, half of the curvature (scalar) power spectrum 
at the end of inflation is due to the sourcing of curvature 
modes by isocurvature modes. 



but this is not a strict requirement per se. However, 
if the total sourcing is tiny (sin A at <C 1) or is very large 
(cos A^r <C 1), then the trigonometric terms will usually 
prevent the bound in equation (|59[) from being satisfied. 
We add that to achieve moderate sourcing, the larger Tss 
is during inflation, the smaller the turn rate must be, and 
vice versa. For reference, Figure 1 shows the value of the 
trigonometric factor sin 2 A jy cos A 2 N as a function of the 
total mode sourcing, Tns- Its maximum value is 0.25, 
which occurs at T-rs = 1- 

Although we have not shown it here and the calculation 
is more difficult, similar qualitative conditions hold for 
general multi-field inflation; the main difference is that 
for multi-field inflation, the sourcing term analogous to 
Tus is a vector, rather than a scalar. This explains why 
it has been difficult to find multi-field models of infla- 
tion that produce large non-Gaussianity: it is not easy 
to find inflationary scenarios that are so dramatically sen- 
sitive to the initial conditions and that involve moderate 
sourcing. The interesting corollary of this is that some 
degree of fine-tuning is needed. Fine-tuning is needed 
both to produce a potential where the mode sourcing is 
so sensitive to the initial conditions and to start in the 
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e-folds of inflation, N 

FIG. 2: For |/jvl| to be large, the amount of sourcing of curvature modes by isocurvature modes (T-jzs) must be extremely 
sensitive to changes in the initial conditions and usually the amount of sourcing must also be moderate. Above is an example 

of a trajectory (solid lines) that meets these two criteria. We use the potential V((/>i,(j>2) = |e _A ^ 2 m 2 </>i, which Byrnes 
et. al. thoroughly investigated in [4q[ . We set A = 0.05 and illustrate the results for two trajectories: (1) one that starts at 
(</>*> <f>2) = (17, 10 -4 ), follows along the ridge, and turns ever so slightly at the end of inflation (solid lines), and (2) a neighboring 
trajectory that starts only |A0J = 0.01 away in field space, but that eventually rolls off the narrow ridge (dashed lines). The 
plot on the left shows the inflationary potential as a function of the fields, along with the two inflaton trajectories. The plots 
on the right show the turn rate, the relative amplitude of isocurvature modes (Tss), and the total amount of sourcing (Tjis) 
as a function of N. Here, the approximately 50-fold difference in the total sourcing stems more from the difference in the turn 
rates for the two trajectories, which both possess large isocurvature modes during all of inflation. Interestingly, the trajectory 
that rolls along the ridge (solid lines) produces |/jvx| ~ 10 2 , while the neighboring trajectory (dashed lines) corresponds to 
\/nl \ ~ 1 [4g|, visually illustrating the role of fine-tuning in achieving large non-Gaussianity. 



very narrow subset of initial conditions where both of the 
above conditions are satisfied. 

In Figure 2, we provide an example of an inflationary 
scenario that meets the two criteria associated with the 
bound in equation (|59[) . In the example below, an unsta- 
ble ridge in the potential provides the perfect conditions 
for making the total amount of sourcing (T-j^s) so sensi- 
tive to the initial conditions. The trajectory of interest 
(solid lines) rolls along the ridge and turns ever so slightly 
at the end of inflation, resulting in moderate sourcing and 
hence |/jvl| ~ 10 2 - By comparison, a neighboring tra- 
jectory (dashed lines) rolls off the ridge, experiences ex- 
tremely strong sourcing, and produces |/jvl| ~ 1, which 
is just below the detection threshold for CMB experi- 



ments. 



B. Calculation of el V*Ttj 5 

Now we proceed to take the calculation of e^j_ • V*Ttj5 
as far as possible. The work we present in this section is 
most applicable to analytically solvable models. We also 
note that in this section, we explicitly indicate when the 
transfer function Tss and certain other quantities are to 
be evaluated. In particular, we use the superscript e to 
denote that a quantity is to be evaluated at the end of 
inflation. 

Calculating the term • V^T-^s is more difficult be- 
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cause Tns is often not a conservative function and be- 
cause the field values at the end of inflation depend on 
the held values at horizon exit. To act the operator V* 
on the expression for Tns hr equation (1321) . we change 
variables and rewrite the equation as a line integral ex- 
pression of the fields. We can convert equation (|32|) into 
a line integral by working in the SRST limit and by re- 
placing the function a = 2^- with its SRST counterpart 
— 2A/||j_, which yields 

Tns = -2 [*' Tss{^4>) (60) 
J<t>„ v2e 



If the integrand of Tns is the gradient of a function, 
then operating V* on Tns simply returns the integrand 
evaluated both at horizon exit and at the end of infla- 
tion, with the latter being times a matrix representing 
the sensitivity of the final field values to the initial field 
values (by virtue of the Chain Rule of calculus). How- 
ever, in general, the integrand of Tns will not be the 
gradient of a function. To account for this, we introduce 
a model-dependent function 7 to represent how much the 
integrand of Tns i n equation (|60p deviates from being the 
gradient of a function. Now operating V* on Tns an d 
using the new function 7, we obtain 



V.Tns = 



2Me 



7 e J 



X 



2Me 



-7 e -L 



T§ s +T ns V* (In 71s), 



<p e — const ' 



(61) 



r 



where 



X> 



dC 



dC 



(62) 



The matrix X arises in the above expression due to the 
Chain Rule of calculus. It captures how a change in the 
initial conditions at horizon exit affects the final values of 
the fields at the end of inflation. Since only changes in the 
initial field vector that are off the trajectory will affect 
the final field values, the variable C parametrizes motion 
orthogonal to the given trajectory; in other words, C is 
constant along every unique trajectory of motion. The 
last new quantity we introduce in equation (|6ip is the 
term V* (lnTJ s )^ =const > which means to take the gra- 
dient of In TJ 5 while holding the amplitude of the isocur- 
vature modes at the end of inflation constant. This term 
can be thought of as some sort of measure of the sensi- 
tivity of TJ 5 to the initial conditions. This term arises 
from the fact that acting V* on Tns in equation (f60|) 
involves differentiating under the integral, which is nec- 
essary since Tss{4>*, <P) depends on 

Considering equation (IBTj) . the e| component is a 
model-independent expression simply because it can be 
related to the time-derivative of Tns > which has a model- 
independent form. In fact, we already included this term 
in the first term on the right-hand side of equation (|58|) 
for Jnl, and hence we do not need to consider it fur- 
ther. Thus, we only need to consider the compo- 
nent of V*Tns- We emphasize that the form we as- 
sume for this component is best applicable to analyti- 
cally solvable models and models where the coupling term 
I Q^yfoy I C 1 or is approximately constant, as will be- 
come clearer later. For non-analytic models, • 'V^Tns 
can be evaluated numerically either directly from the ex- 
pression for Tns ° r y i a an alternative expression that we 
present later in this section. In the remainder of this sec- 
tion, we discuss the three terms X, TJ 5 , and 7, which 



arise in the expression for ■'V^Tns- We first show that 
X has a model-independent form. Then, we discuss the 
transfer function TJ 5 , how it affects Tns an d the gradi- 
ent of Tns 1 an d how it determines the model-dependent 
term 7. Along the way, we consider the exact solutions 
for product and sum potentials, casting previous results 
for /jv l for these models 0, |45| in a more geometrically 
and physically transparent form. 

Now we find the matrix X and prove that it has a 
model-independent form. First, since C is constant along 
a given trajectory, 



C = 4! ■ VC = 0. 



(63) 



Therefore, VC must be proportional to ej_ — that is, 



VC = |VC|e_L. 



(64) 



Next, consider -jjjr- Since C is constant along a trajec- 
tory, a change in C corresponds to motion orthogonal 
to the trajectory, and hence must be parallel to e e ± . 
Now combining this fact with equation (1641) and with 



dC 
dC 



d<{> e 
~dC 



V„C 



(65) 



and equation ([M)) into equation (|6"2")l yields 



implies that t5t = |VC|j 1 e5_. Substituting this result 



(66) 



X = 



ivci 

IVCL 



.(el) 3 



Finally, we express the ratio of the norms of the gradients 
of C at horizon exit and the end of inflation in terms of a 
physical quantity: the relative amplitude of the entropy 
modes. Since SC = 5<fi ■ VC, then for a given variation 
in the trajectory, SC, we have 



IVCI 



= IVCL 



(67) 
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Combining equations (|66[) and ([67 
the model-independent expression 



we finally arrive at 



^T| 5 el(ei) 2 



(68) 



This interesting result shows that the sensitivity of 
the final field values to the initial field values can be 
given very simply in terms of the relative growth or de- 
cay of entropy modes. In other words, the evolution of 
entropy modes mirrors whether neighboring trajectories 
converge or diverge over time. In scenarios where neigh- 
boring trajectories converge ("attractor solutions"), the 
entropy modes decay. However, when neighboring trajec- 
tories diverge, the entropy modes grow. That such a rela- 
tionship should hold between the convergence / divergence 
of neighboring trajectories and the evolution of entropy 
modes makes sense. From a geometrical perspective, we 



intuitively expect that a positive curvature along the en- 
tropic direction focuses neighboring trajectories, whereas 
a negative curvature creates a hill or ridge in the po- 
tential, causing neighboring trajectories to diverge. But 
we also know that the curvature along the entropic di- 
rection determines the evolution of entropy modes. By 
equation (|28[) . the entropy modes grow when M±± < 
and decay when M_l_l > 0, and how quickly they do so 
depends on the magnitude of the curvature. Combin- 
ing these two facts together, we could have concluded 
that the divergence/convergence of neighboring trajec- 
tories must correlate with the growth/decay of entropy 
modes, without even deriving this result. Nonetheless, 
equation (|6"8)l gives the precise relationship explicitly. 



Now substituting equation (|68l) into equation (l6Tj) and 
projecting the result onto ej_, we obtain 



J 



V,T K s = (2M 



_L_L 



/2iT 7 *) - (2MLl + V^Tele) (Tj 5 ) 2 + V2T*T ns e* ± ■ V, (In T$ s ) =c 



(69) 



T. 



sax v e cos ( 



55 



The above equation shows that the sensitivity of Tns to for these models can be approximated by 
the initial conditions is determined by M±±, e, TJ S , Tns, 
and the model-dependent factor 7. The above relation 
gives us a nice way to understand when the sourcing is 
very sensitive to the initial conditions, based on the geo- 
metrical and physical attributes of an inflationary model. 

Finally, we consider the model-dependent quantity 7. 
Recall that we defined 7 so that it is zero whenever the 
integrand of Tns is the gradient of a function. This oc- 
curs for product potentials, defined as 



(73) 

sin ft* cos 

Substituting = Q 1 anc } equation (|73|) into equation 
© yields 



Tns = 



26' IT 006 ' dN, 



N, 



V = V 1 (<h.)V2(<h) ) 



(70) 



1 



sin V* cos ( 



and can be attributed to the fact that in these models, 
the two fields evolve independently of each other. 7 is 
therefore non-zero whenever the evolutions of the two 
fields influence each other. 

We can see that 7 is zero for product potentials as 
follows. For product potentials, the isocurvature mass 
equals 



smtf* cost/* J Nt dN 
which integrates to give 



(sin 2 6) dN, 



(74) 



tan0* + tan e TJ s 



sin tf» cos ( 



sin 2 6* 



(75) 



M, 



llll 



M, 



(tan0-cot0)M||_L, (71) 



Now we take the gradient of the above transfer function 
and use that 



where is the polar coordinate in the 



plane, 



tan0 



(72) 



Ve|| = V(cos0,sin0) = _ M f^ e T 



'2c 



(76) 



Using —Mim ~ ^ = 0' an d plugging equation (TfTT) into 
equation (|32p . one finds that the transfer function TJ 5 



in the SRST limit for any two-field model of inflation. 
Finally, projecting the result onto e* ± and using equation 
flTTl) . we find 
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/2i7el • V^Tns = 2M* ±± - 2M e ±± (T§ s ) 2 + V^T KS e* ± • V* QnTJ 5 )^ 
= [2 + (cotfl, - tan^)T K5 ] M xx - 2M xx (Tf s ) 2 , 



M* 



A/ 



ll-L 



M xx - 2M xx (Tj 



55 ; 



(77) 



where again TJ S for product potentials is given by equa- 
tion (|73]l. Comparing equation (|77| to equation (joTT) in- 
deed shows that 7 = 0. 

Now for the general case of two-field inflation, we can 
use a similar procedure to find the model-dependent term 
7. As the transfer function Tss determines how much the 
integrand of T-rs deviates from being the gradient of a 
function and hence determines 7, we start by considering 
Tss- We begin by finding a general expression for the 
transfer function Tss- Starting from equation (42) in 
[85j |. after some algebra, we can show that this implies 



that TJ 5 takes the form 



55 



sin 9 e cos 9 e 
sin 6* COS0* 

sin 9 e cos 9 e V e 
^sin 9* cos 9*V* 

where A/ 12 = 8182 In V and 
s(0„,0 e ) = exp 



exp 



jv, 



-dN 



N, 



s(0„0 e ) 



y<9i<9 2 y 



(78) 



V In V • d0 



(79) 



This means that whenever ( g^/^y ) ^ s a cons tant, 
s((p lf , </> e ) becomes analytic and hence Tss becomes ana- 
lytic. For product potentials, M12 = 0, reproducing the 
result we derived in equation (|73l) . For sum potentials, 
defined as 



V = V 1 (cb 1 ) + V 2 (<f )2 ), 
the coupling term 8182V — 0, and so 
sin 9 e cos 9 e V e 



rjie 

J 55 — 



sin 9* cos 6>*y» 



(80) 



(81) 



Equation ([8~Tj) for TJ^ can also be used as an approxi- 

•C 1 during all 



L ss 

mation for scenarios in which 



Vdid 2 V 
d!Vd 2 V 



of inflation. In the more general case where 



diVd 2 V 



approximately constant during inflation, TJ 5 can be ap- 
proximated by an analytic function similar to equation 
(IBT1) . but possessing additional powers of (V e /V r ). And 

whenever ^^75^7^ is approximately constant — which 

is automatically true for all product and sum potentials 
— the following relation holds: 



/2eT el -V* (In T| s )0 =const = (cot0» -tan0*)M 



_U_! 

(82) 



as gradients of functions of V do not contribute to 
e x • V*(lnTJ <s )0 =CO nst- Otherwise, when the term 

( "diV% v ) ^ s n0 ^ a cons tant, this adds extra terms to 
\/2e7e ± • V*(lnT| s )^ e=const that contribute to f NL . 

Now plugging the general expression for Tss into equa- 
tion (|32p and integrating by parts again using the fact 
that (sin 2 9)' = 2 sin cos 00', we obtain 



T 



RS 



tan 0* + tan e Tj 5 



1 



sin 0, cos 0*y 



* ' * JJV, 



(83) 

sin 2 0-£- (Vs(N*,N))dN. 



dN 



Now the perpendicular component of the gradient of T-rs 
can be calculated directly using the above equation, 4 
where equation (|76[) comes in handy. In particular, for 
sum potentials, the integral in equation (|83[) evaluates to 
yf - y 2 * , yielding 



Tns = - tan ( 



Vo* 



sin 9* cos0*y* 



tan 9 e 



\7 



sin 9 e cos 9 e V e 



n 



55- 



(85) 



For all other potentials, the integral in equation (|53")) typ- 
ically must be computed numerically, which is still usu- 
ally the fastest route to finding e* • V*T-r.s. However 
we instead use the above expression to find the model- 
dependent term 7. Taking the gradient of equation (|83[ 
and factoring it into the form of equation (1611) , we con- 
clude that 7 can be written as 
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4 From equation 1 183 II , we can also derive an upper limit for Tjiq 
for general two-field inflation whenever the SRST limit is a valid 
approximation. Using the fact that sin 2 9 < 1, we find that from 
equation (183 P . 

T-rs < cot 6* - cot e T§ s . (84) 



7e 



tan 8 e e e ± ■ V e ln(s(0 t , e )) 



1 



sin 0* cos d*V* 



sin 2 V(l/ S (</>„, </>))-g!0 



(86) 



- 0^— constant 



where the gradients are taken with respect to holding the 
initial field vector constant. The resulting expression for 
7 at horizon exit has the same functional form, of course. 

For sum potentials, the function s(4> !¥ ,<p) vanishes, 
leaving only the perpendicular component of the gradi- 



ent of the integral JJ e sin 2 9 V V ■ d<j> = V£ - V{ . Hence 

we find that 7 = ~\[2e. Therefore, the expression for 
e*j_ ■ V*T-r,s for sum potentials can be written as 



/^el • V*T ns =(2M* ±± - 2e„) - (2M e ±± - 2e e )(TJ 5 ) 2 + V^T ns e* ± ■ V* (lnTj s ) 0e=c 
= [2 + (cot0» - tanflOTrcs] M* ±± - 2e* - (2M e ±± - 2e e )(T^) 2 , 



M l±_ ~ M \\\\ - 2e * 



T. 



11s 



Ml 



2e,-(2Afi ± -2e e )(Tj iS ) 2 , 



(87) 



where TJ 5 is given by equation (|8ip. 

For all other potentials, equations f|69[) . ([79|l , and (|86|) 
represent the prescription for calculating e* ± ■ V^T^g. 
For weak coupling among the fields, we expect 7 to be of 
order the slow-roll parameters. However, this term may 
be larger in the limit of strong coupling. 



C. Conditions for Large /jvl 

As we showed in Section pV A|) . if the power spec- 
tra are nearly scale- invariant, the magnitude of /jvl can 
be greater than unity only if | sin A 2 ^ cos Afjy/2e* e* ± ■ 
V *Tjis\ > 1. Satisfying this bound requires that two 
conditions be met: (1) that T-rs be extremely sensitive 
to changes in the initial conditions perpendicular to the 
given trajectory and (2) that the amount of sourcing 
be non-zero. Typically, it also means that the sourc- 
ing of curvature modes by isocurvature modes must be 
moderate, although this is not strictly required; rather, 
very weak sourcing (sin Ajv < 1) or very strong sourcing 
(cos A^r 3> 1) simply makes the bound in equation (|59|) 
extremely hard to satisfy. 

The first condition can be understood very simply: in 
order to produce a large degree of skew in the primor- 
dial fluctuations, perturbations off the classical trajec- 



tory must move the inflaton onto neighboring trajecto- 
ries that experience very different dynamics for the field 
perturbations. To satisfy this condition requires that the 
profiles of the turn rate and/or of the relative amplitude 
of isocurvature modes (T55) be dramatically different for 
neighboring trajectories. The second condition, that the 
sourcing must be non-zero, requires that the turn rate 
not be zero for all of inflation. And to achieve moderate 
sourcing, which is usually needed to satisfy the bound in 
equation (|59p. the larger T$s is, the smaller the turn rate 
must be, and vice versa. 

One question that naturally arises is whether combin- 
ing the need for moderate sourcing with the requirement 
that the total sourcing be dramatically different among 
neighboring trajectories gives us any constraints on the 
possible ways that the turn rate and/or T$s can vary and 
still produce large |/jvl|- The answer to this question is 
yes. We can best see this by considering the sourcing 
function T-ks- Consider the case where the amplitude 
of isocurvature modes never exceeds its value at horizon 
exit, i.e., Tss{N*,N) < 1. In this case, T-rs can never 
exceed 

r ns <2 e'dN = 2(e e -e*), (88) 

Jn, 

where again is the polar angle for the field velocity 
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vector. For the most common scenarios, the field ve- 
locity vector does not turn through an angle of more 
than 90°, yielding a bound of Tjzs < tt. Let us com- 
pare this bound of Tns < tt with a numerical example. 
If e* = 0.02 and we assume nearly scale-invariant scalar 
and tensor spectra, then we need e^j_ • V*T-r.s > 60 in 
order to produce \Jnl\ ~ 3. This seems extremely dif- 
ficult to achieve given the bound of Tjis < tt and the 
need for moderate sourcing. Now if within this set of 
scenarios, we consider the subset where M±± > dur- 
ing all of inflation, then neighboring trajectories must 
converge (or at least not diverge) over time. Since neigh- 
boring trajectories remain close to each other during all 
of inflation, Tns cannot differ widely among neighboring 
trajectories without discontinuous or other extreme fea- 
tures in the potential that violate the SRST conditions. 
To be more precise, this would require the speed-up rate, 
turn rate, and/or entropy mass to be hugely varying in 
the direction orthogonal to the given trajectory, which 
effectively constitutes a violation of higher-order SRST 
parameters. And on an intuitive level, it is not possi- 
ble for a neighboring trajectory to have a much larger or 
smaller turn rate for a substantial period of time with- 
out having the two trajectories diverge, 5 which violates 
M±± > 0. Therefore, since large non-Gaussianity can- 
not be produced under these conditions, either Tss > 1 
and/or Mj_j_ < at least sometime during inflation. 

In Section (jIVBp . we took the calculation of V2e* • 
V»Tk5 as far as possible, trying to better understand 
when this term is large in magnitude in analytically solv- 
able and similar models. From equation (1691) . we found 
that y/2e* e* ± ■ ^ *Tjis will be large in magnitude if at 
least one of the following three conditions is met: 

5 We make the usual unstated assumption that there are no clas- 
sical degeneracies in the gradient of In V, which means that tra- 
jectories cannot cross each other. 



1. 2M^_ ± — ^2e»7» is large in magnitude, 

2. (2M^_ ± — -\/2e e 7 e )(TJ 5 ) 2 is large in magnitude, or 

3. T| 5 is very sensitive to changes in the initial con- 
ditions orthogonal to the inflaton trajectory. 

In conventional slow-roll, the magnitudes of Mj_ ± and 
the other slow-roll parameters at horizon exit are sig- 
nificantly less than unity, so only the latter two condi- 
tions can be satisfied. If we assume conventional slow- 
roll at horizon exit and additionally that the magnitudes 
of MJj_ and 7 e are at most of order O(l)-O(10), then 
the second condition above requires that TJ 5 is at least 
of order unity. If we tighten the constraints even further, 
requiring that |Af^jJ <C 1 during all of inflation, then the 
second condition becomes even more stringent, requiring 
that Tg S be very large, at least of order O(100). 

Interestingly, we can show that essentially the same 
conditions for large \/2e*e* L ■ V»T-r.s arise via an alterna- 
tive approach. Anytime VJV^ ^constant — F — which 6 
includes all scenarios in which N is a function of only the 
initial and final fields — we have 

V*iV = -F» + ATF e , (89) 

and hence by equations ([4*5]) and ([55]) . 



Tns — - 



2eF, 



V2eFA TJ 5 



(90) 



Now taking the perpendicular component of the gradient 
of Tjis, we obtain 

6 By this expression, we mean that the gradient of N evaluated 
when holding the final field vector constant equals F. 



J 



'2e* V*i Tns — — 



2eV 



(VTeF±) 



2eV 



2eF, 



r 



2eFj_ 



(91) 



where Vj_ = ej_ ■ V. So under similar assumptions, it 
again appears that the amplitude of isocurvature modes 
at the end of inflation and the sensitively of TJ^ to a 
change in initial conditions orthogonal to the given clas- 
sical trajectory largely control the magnitude of /jvl- 

Let us now consider the above conditions in the context 
of the two most popular categories of models: product 
potentials and sum potentials. Expressions for /jvl for 
product potentials and sum potentials were first found by 
HU and [HI, respectively. Our versions of the same ex- 
pressions — equation ([58]) coupled either equation ([77]) 
or (|87l) . respectively — give somewhat more direct in- 
sight into how Jnl depends on the physical and geomet- 
rical features of an inflationary model. The conditions 



for large bispectra in both product and sum potentials 
were discovered by Byrnes et. al. 46]. For product 
potentials, they found that /jvl will be large in mag- 
nitude when either one of the two fields starts with far 
more kinetic energy than the other (either cot ^> 1 or 
tan 6* 3> 1) and when the asymmetry in the kinetic ener- 
gies of the two fields diminishes significantly by the end 
of inflation. The reason that these two conditions pro- 
duce large bispectra in product potentials is that they 
together guarantee three things: that TJ 5 ^ 1 (which 
follows from equation ([75]) ). that TJ 5 is very sensitive to 
changes in the initial conditions orthogonal to the infla- 
ton trajectory, and that the turn rate is small yet signifi- 
cant enough to produce moderate sourcing. Interestingly, 
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this means that for product potentials, both the second 
and third conditions for large -y/2e» • V*T-r,s are always 
simultaneously satisfied, as the third condition combined 
with the requirement that the total sourcing is moder- 
ate guarantees the second condition, and vice versa. For 
sum potentials, the conditions for large bispectra are a 
bit more complicated to untangle but end up being sim- 
ilar; however, the expression for TJ S contains a factor 
of t£-. Although we state the conditions for large \Jnl\ 
differently and slightly extend them in scope by relax- 
ing constraints on M±±, these conditions otherwise agree 
with those uncovered by Byrnes et. al. [46| . 

As we concluded in Section pV A[) that large non- 
Gaussianity requires the profiles of the turn rate and/or 
T$s to be dramatically different for neighboring trajecto- 
ries, we might worry why we did not uncover any explicit 
conditions that involve the turn rate. There are a few rea- 
sons for this. First, the lowest-order time variation of the 
turn rate is a model-dependent function of the isocurva- 
ture mass, ft = Miiii — M±± and of the coupling between 
the fields, so the turn rate is implicitly included in the 
above conditions. Second, this should not worry us as the 
turn rate is constrained not to be too large by the con- 
straints on scale- invariance and is constrained not to be 
too small by the need for moderate sourcing. And third, 
the difference in the turn rate between two neighboring 
trajectories cannot be large for a sustained amount of 
time without neighboring trajectories diverging, so M±± 
itself also tells us whether larger variations in the turn 
rate between neighboring trajectories are possible. 

Now we consider the geometric implications for the in- 
flationary potential. We argued earlier that it is not pos- 
sible to achieve large non-Gaussianity during inflation if 
both T S s(N*,N) < 1 and M±± > during all of infla- 
tion. Therefore, one or both of the conditions must be 
violated sometime during inflation to produce large non- 
Gaussianity. If we make the conventional assumptions 
that e* <tC 1, e increases significantly (but not necessarily 
monotically) in order to end inflation, and that e never 
drops below its value at horizon exit, then the only way 
to satisfy Tss{N*,N) > 1 is for the entropy modes to 
grow at some point during inflation. 7 As we showed ear- 



If we relax the assumption that e never decreases below its value 
at horizon exit, then it may be possible to produce large |/jvlI 
in two-field inflation without any negative curvature along the 
entropic direction. Indeed, Byrnes et. al. |61| showed that this 
is possible for a three-field model: a two-component hybrid in- 
flation model, where we are counting the waterfall field as the 
third field. During the first phase of inflation, which is governed 
by a vacuum-dominated sum potential and during which only 
two of the fields are active, e exponentially decays after modes 
exit the horizon. The exponential decay of e causes the isocur- 
vature modes to grow dramatically during the initial phase of 
inflation, which produces large non-Gaussianity even before the 
waterfall field comes into play. This large exponential decay in e 
then circumvents the absolute need for a negative entropy mass 
at some point during inflation. However, for this model to then 



lier, the entropic modes grow when the curvature along 
the entropic direction is negative, or equivalently, when 
neighboring trajectories diverge. Hence under conven- 
tional assumptions, the two conditions become one and 
the same, and we require that M±± < at least some 
time during inflation. Geometrically, this means that the 
inflaton must roll along a ridge in the potential for some- 
time during inflation — that is, the inflaton trajectory 
must be unstable. Perturbations orthogonal to the clas- 
sical trajectory must result in neighboring trajectories 
that diverge, giving rise to widely different inflationary 
dynamics. Figure 2 illustrates one geometric realization 
of a potential that meets these criteria. Moreover, the 
potential must not only possess a ridge, but the initial 
conditions must be fine-tuned so that the inflaton rolls 
along the ridge for a sufficiently long time. Conversely, 
whenever attractor solutions exist, non-Gaussianity will 
be small. 

Also in the geometric picture, it is necessary that the 
inflaton trajectory turn somewhat during inflation, as we 
argued earlier that \Jnl\ can only be large if the total 
sourcing of curvature modes by isocurvature modes is 
non-zero and usually it must also be moderate. Other- 
wise, if the sourcing is too weak or too strong, then the 
trigonometric functions representing the sourcing effects 
will still force |/atl| to be small. Interestingly, how much 
of a turn is needed depends on the relative amplitude of 
the isocurvature modes during the turn. If the isocur- 
vature modes are small at the time, then a larger turn 
is need. However, if the isocurvature modes are large, 
then only a minuscule turn in the trajectory is needed to 
produce moderate sourcing. 

Cumulatively, these results explain why it has been so 
difficult to realize large non-Gaussianity in two-field in- 
flation. The potential must be fine-tuned enough to pos- 
sess a steep ridge, while the initial conditions must be 
fine-tuned enough that the inflaton rolls along the ridge 
for a significant length of time during inflation, but also 
slightly turns. Nonetheless, a few two-field scenarios that 
produce large non-Gaussianity have been identified. Hy- 
brid inflation, where the waterfall field counts as one of 
the two fields, is a great example of an inflationary sce- 
nario possesses a negative entropy mass during inflation 
and that meets these conditions. As a second example, 
Byrnes et. al. [46[ showed that product potentials where 
one of the two fields dominates the inflationary dynam- 
ics (meaning that the kinetic energy of one field is much 
larger than the kinetic energy of the other) during all of 
inflation, but the subdominant field picks up speed log- 
arithmically faster than the dominant field also produce 



to be viable, a third field, the waterfall field with its associated 
negative mass, is needed to end inflation. We also note that 
double inflation models with very high mass ratios do violate the 
assumption that e never drops below e* and hence exceed the 
bound Tfzs < TTj but nonetheless, they do not to produce large 
non-Gaussianity. 
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large |/jvl|- Finally, we note that in investigating sum 
potentials, many have found it more difficult to find sce- 
narios in which |/jvx| is large, but the primary reason 
for this trouble is that the focus has been on scenarios 
in which the entropy mass is strictly positive. Therefore, 
observable large non-Gaussianity should be more readily 
achieved in sum potentials by searching for those scenar- 
ios that allow the entropy mass to be negative. 



V. TRISPECTRUM 

Now we calculate the local trispectrum. The local 
trispectrum can be expressed in terms of two dimension- 
less non-linear parameters, tnl and gNL'- 

T n = t nl [TMIki + k 3 \)V n (k 3 )VTz(ki) + 11 perms] 
54 

+ ^g NL [Pn{k 2 )Vnik 3 )V K (ki) + 3 perms] . (92) 
where under the SN formalism [86l [8 



servables, we focus on tnl- Equation (J93J) for tnl is 
equivalent to 



tnl 



|ViV| 2 



(94) 



Using equation (IM|) and the fact that we already calcu- 
lated the form of the vector v ^^| 2 ew in equation (|54| , 
we can quickly arrive at the answer. Dividing equa- 
tion (I56[) by cosA^v, we obtain the eii component of 

VV T jVe N . 



|VW| 2 



eiiVV^iVeAr 1 



|V7V|2 



= - cos A N (n-ji - n T ). (95) 



Similarly, dividing equation (|5T[) by sulAat gives us the 
ej_ component: 



tnl 



54 

25 9NL 



e N VV T NVV T Ne N 
' |ViV| 4 ' 

e N e N eJj V VV7V 



|VJV| 2 



(93) 



e ± W 1 Ne N 
\VN\ 2 



= COS A /v X 



(96) 



M: 



l-L 



sin A n cos AnV%£* e ± ■ V*T K5 



Since the expression for gNL is less illuminating and 
cannot be completely expressed in terms of other ob- 



Substituting the above two equations into equation (|94|) , 
we find that 



J 



tnl = | cos 2 A N {nn - n T ) 2 + cos 2 A N 



M||_l + sin A N cos AnV^^I ■ VT RS 
I 



(97) 



Using equations (|37|) . (J40J) , and (|58|) . we can write tnl consistency condition: 
completely in terms of observables, giving the following 



tnl 



(l - r 2 c ) [n n - n T } 2 + 



1 

ZT 

1 c 



6 1 

-Inl + ^ (! - r c) ( n n ~ n T ) 



(98) 



where recall that tq is the curvature-isocurvature corre- 
lation. This gives us a new consistency relation that is 
unique to two-field inflation and that relates the observ- 
ables tnl, Snl, n n , n T , and r c . 



was shown [671 ] that 



TNL > 



tInl 
5 



(99) 



Examining equation (|98[) , since the magnitudes of n-ji 
and tit are constrained to be much smaller than unity, 

tnl can only be large if f'^r) is large. Previously, it 



proving that tnl will be large whenever f NL is large, but 
this result still left open the question of whether it is pos- 
sible for t N l to be large if ff iL is not. This question has 
since been answered affirmatively for particular models 
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(e.g., 0, [6a|). Above we show that it is more generally 
possible for tnl to be large even if f NL is not, but only if 
r c is small. As rc = sin A jv reflects the degree of sourc- 
ing of curvature modes by isocurvature modes, this shows 
that the commonality in two-field models where tjvl is 
detectably large when /jvl is not is that the sourcing 
effects (the multi-field effects) are weak. 

In the limit where |/jvl| ^ 1, equation ([95]) reduces to 



tnl 



G 



I JVL 



(100) 



In this limit, how much tnl exceeds the Suyama- 
Yamaguchi bound in equation (|99[) depends only on rc- 
While we might naively expect that making rc as small 
as possible would maximize the value of tjvl, this is not 
necessarily the case. This is because there is a trade-off: 
scenarios in which the multi-field effects are very small 
{rc <C 1) behave in many ways like single-field mod- 
els and hence they are likely to produce small values for 
| /jvl |- Therefore, reducing rc even further while preserv- 
ing larger values for |/jvl| typically comes at the expense 
of even more fine-turning. Therefore, for detectable non- 
Gaussianity without excessive fine-tuning, we might ex- 
pect more typically that tjv l is no more than one to two 
orders of magnitude greater than f NL . 

Using a similar approach, an expression for gNL can 
be found. However, the result can be expressed only 
partially in terms of observables. Below, we cast the 
result in the most simple and transparent way. Starting 
by operating r^m e N ■ V on the expression for /jvl, and 
then using the definitions of the non-linear parameters 
and that 



tnl 



e T N V(\VN\e T N )VV T Ne 



N 



|VA| 4 

e^V(|VA^|)e^ VV T Ae,v 



|VA| 4 
e N (Ve N )VV T Ne N 
IVAl 3 



(101) 



gNL can be written as 



54 „ ,(§r 

7j^9NL = -2t nl +41 -/jvl 



/ r T ( 6 

T ejv-V(--/jvL 



(102) 



Hence gNL can be large in magnitude only if tnl is large, 
f NL is large, and/or if /jvl varies dramatically in a small 
neighborhood about the initial conditions. Therefore, if 
gNL is large, but neither tnl and /jvl are, then it means 
that the inflaton trajectory is near neighboring trajecto- 
ries that do produce large /jvl and/or that /jvl has very 
strong scale-dependence. 



VI. CONCLUSIONS 

In this paper, we have derived formulae for the local 
form of the bispectrum and trispectrum in general two- 
field inflation. In particular, we found semi-analytic ex- 
pressions for /jvl, tnl, and gNL, the only non-linear 
parameters in the expressions for the bispectrum and 
trispectrum whose magnitudes have the potential to be 
at least of order unity. To do so, we worked within the 5N 
formalism, which expresses the bispectrum and trispec- 
trum in terms of gradients of N, where N is the number 
of e-folds of inflation. To perform the calculation, we in- 
voked the slow-roll and slow-turn approximations, and we 
used a unified kinematical framework, the transfer ma- 
trix formalism, and a general expression for the evolution 
of isocurvature modes. 

We showed that /jvl can be written in terms of sines 
and cosines (related to the degree of sourcing) times 
n-ji = 1 — n s , riT, the turn rate at horizon exit, and 
\f2ele* L ■ V*T-r,s, where T-rs is the transfer function 
that encodes the relative degree of sourcing of curvature 
modes by isocurvature modes. As the magnitudes of all 
quantities but the term sin 2 Ajv cos 2 Ajv-\/2e* • V*T-r,s 
are constrained to be less than unity, |/jvl| can only be 
large when (1) Tus is extremely sensitive to a change 
in initial conditions orthogonal to the inflaton trajectory 
and (2) the total sourcing is non-zero, though usually 
the total sourcing must also be moderate. The former 
condition makes sense on an intuitive level, as to pro- 
duce a large amount of skew in the primordial pertur- 
bations, fluctuations off the classical inflaton trajectory 
must result in very different inflationary dynamics for 
the field perturbations. Now since Tns is an integral of 
the turn rate and the relative amplitude of isocurvature 
modes (Tss), the former condition implies that neigh- 
boring trajectories must have dramatically different turn 
rate profiles, T ss profiles, or both. Though we only pre- 
sented proofs of these conditions for two-field inflation, 
similar conditions hold for multi-field inflation as well. 

Next, we found an expression for ^/2e^e* ± ■ V*T-r.s for 
analytically solvable and similar scenarios that depends 
on the entropy mass, the isocurvature transfer function 
TJ 5 , and a model-dependent correction 7, which quan- 
tifies the coupling between the fields. Invoking mini- 
mal assumptions about the terms in equation (|6T))) , we 
showed that for y/2e*e* ± ■ V *T-ns to be large requires 
that Mxj_ — 7 be large at horizon exit, that the rela- 
tive amplitude of isocurvature modes (Tss) at the end 
of inflation be large, and/or that T| 5 be very sensitive 
to changes in the initial conditions perpendicular to the 
inflaton trajectory. 

We then further explored the conditions for large non- 
Gaussianity in general two-field inflation. After proving 
an upper bound for Tus in the case where Tss < 1 
during all of inflation, we argued that if neighboring 
trajectories do not diverge, then due to constraints on 
higher-order SRST parameters, the amount of sourcing 
cannot vary dramatically among neighboring trajectories 
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and hence non-Gaussianity cannot be large. Therefore, 
either T55 > 1 or Mj_j_ < sometime during inflation, 
and under some minimal assumptions, these conditions 
become one and the same. Geometrically, this means 
that Jnl will be large only if the inflaton traverses along 
a ridge in the inflationary potential at some point dur- 
ing inflation and the inflaton trajectory at least slightly 
turns so that the total sourcing of curvature modes by 
isocurvature modes is moderate. Unfortunately, though, 
this implies that some fine-tuning of the potential and/or 
the initial conditions is needed both to produce a steep 
enough ridge and/or to situate the inflaton on top of 
the ridge without it falling off too quickly and yet still 
slightly turning. Inflationary scenarios that are attrac- 
tor solutions therefore cannot produce large Jnl- This 
explains why it has been so difficult to achieve large 
non-Gaussianity in two-field inflation. Moreover, it ex- 
plains why large non-Gaussianity arises in models such 
as hybrid and multi-brid inflation, axionic A/"-flation, and 
tachyonic (p)reheating. The common denominator of 
these models is a significant negative curvature (mass) 
along the entropic direction. 

Finally, we showed that the calculations of tjvl and 
Qnl are very similar to that of Jnl- tjvl can be written 
entirely in terms of the spectral observables Jnl, n>n, 
tit, and rc, where rc is the dimensionless curvature- 
isocurvature correlation. This provides a new consis- 
tency relation unique to two-field inflation. Moreover, it 
sheds new light on the Suyama-Yamaguchi bound tml > 

(jInl) 2 , showing that for \f NL \ > 1, t nl = ■ 



Though theoretically one could attempt to minimize rc 
to make tml even larger relative to Jnl-, this usually 
comes at the unwanted expense of further fine-turning. 
We also calculated the trispectrum parameter g^L and 
showed that it can only be large in magnitude if t^l is 
large, f^ L is large, and/or Jnl varies dramatically in a 
small neighborhood about the initial conditions. 

Our results for the local bispectrum and trispectrum 
from inflation allow us to better test and constrain 
two-field models of inflation using observational data. 
Our results also provide better guidance for model- 
builders seeking to find inflationary models with large 
non-Gaussianity. In the future, it will be interesting to 
explore the range of shapes of ridges that give rise to 
large non-Gaussianity and to better understand the de- 
gree of fine-turning needed in the potential and/or initial 
conditions. Finally, it is important to better understand 
the model-dependent nature of (p)reheating and the af- 
termath of inflation, to understand the impact on the 
primordial non-Gaussianity from inflation. 
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